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Inhomogeneity in the ground state is an intriguing, emergent phenomenon in magnetism. Recently, 
it has been observed in the magnetostructural channel of the geometrically frustrated Q:-NaMn02, 
for the first time in the absence of active charge degrees of freedom. Here we report an in-depth 
numerical and local-probe experimental study of the isostructural sister compound CuMn 02 that 
emphasizes and provides an explanation for the crucial differences between the two systems. The 
experimentally verified, much more homogeneous, ground state of the stoichiometric CuMn 02 is 
attributed to the reduced magnetoelastic competition between the counteracting magnetic-exchange 
and elastic-energy contributions. The comparison of the two systems additionally highlights the role 
of disorder and allows the understanding of the puzzling phenomenon of phase separation in uniform 
antiferromagnets. 


Although phase separation in a uniform system is a 
widespread phenomenon in diverse fields of matter—^— , 
ranging from biological systems^^— , to soft matter— i^, 
and strongly correlated electron systems^rJ^, in mag¬ 
netism the microscopic pattering has been, until recently, 
almost exclusively limited to thin ferromagnetic (FM) 
fihnsi^di. In this case, such a pattering is a trade-off 
between minimizing the exchange and the dipolar ener¬ 
gies. It thus represents one possible manifestation of a 
general requirement of multiple competing phases that 
can lead to inhomogeneous states. Lately, it has be¬ 
come increasingly apparent that a similar competition 
between energetically nearly equivalent phases is also re¬ 
sponsible for phase separation in geometrically frustrated 
spin systemsi^^— that are generically characterized by 
ground-state degeneracy22. However, the balance be¬ 
tween the competing phases in these systems is generally 
much more delicate and, therefore, poorly understood. 

Recently, the spatially anisotropic triangular antifer- 
romagnet Q;-NaMn02, with dominant intrachain (Ji) 
and geometrically frustrated interchain (J2) antiferro¬ 
magnetic (AFM) exchange interactions (inset in Fig¬ 
ure [Hi), has been highlighted as a paradigm of a phase- 
separated ground state in the absence of active charge 
degrees of freedoroi^. Its AFM order that sets in be¬ 
low the Neel temperature Tn = 45 K is accompanied by 
a simultaneous structural deformation^^. This was ini¬ 
tially suggested as being a phase transition from the high- 
temperature monoclinic {C2jm) to the low-temperature 
triclinic (PI) crystal structure^^. However, more de¬ 
tailed, recent experiments have shown that the magnetic 
order fails to drive this improper ferroelastic transition to 
completions^. Instead, an intricate magnetostructurally 
inhomogeneous state on the nano-scale has been discov¬ 
ered below the Tn- Such a state was suggested to be an 
unforeseen consequence of the subtle interplay between 
the geometrical frustration and the competing structural 


phasesi^. 

In order to fully understand this novel phenomenon, 
further theoretical studies and experimental investi¬ 
gations of related compounds are of paramount im¬ 
portance. In this respect, a comparison with the 
crystallographically^i and magnetically^^ analogous sis¬ 
ter compound CuMn02, known as the mineral crednerite, 
is particularly relevant. Here, in contrast to a-NaMn02, 
the emergent magnetic order below Tn = 65 K is believed 
to lift the macroscopic degeneracy in the spin space com¬ 
pletely, by inducing the monoclinic-to-triclinic structural 
phase transitions^. This spin-induced phase transition 
is witnessed by the splitting of several families of nu¬ 
clear Bragg reflectionsiSiiSl. It was suggested to reflect 
the strong magnetoelastic (ME) coupling that allows for 
the development of shear strain at a low energy costSl. 
Interestingly, the strain is significantly enhancedS^ in 
the off-stoichiometrioS^ Cui+a;Mni_2;02, where Tn is re¬ 
duced and the structural transition temperature is fur¬ 
ther suppressed with increasing x even for small doping 
levels^Si^i^. Moreover, in off-stoichiometric samples the 
interlayer ordering changes from AFM, observed in stoi¬ 
chiometric CuMn02, to FM2^i22i2ij which was attributed 
to a partial substitution of Cu^+ for Mn^+ that should 
effectively change the interlayer exchange coupling from 
AFM to FM22. On the other hand, this implies that 
CuMn02 may be very close to an electronic instability, 
possibly of a similar kind to that found in stoichiometric 
Q;-NaMn02. 

Obtaining in-depth information about the magnetic 
and structural properties of CuMn02 on the local scale 
should clarify the differences with respect to the isostruc¬ 
tural a-NaMn02. Such knowledge would also help 
to address the pending issue of the microscopic ori¬ 
gin of the phase-separation phenomenon in geometri¬ 
cally frustrated magnets. The most obvious ambigui¬ 
ties arise from questions like: why is the Tn enhanced 
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in CuMn 02 compared to a-NaMn 02 , despite the theo¬ 
retically predicted sizably smaller exchange interactions 
in the former compound^; what is the role of the ME 
coupling in establishing the structural distortion below 
the Tn, what is the role of disorder; and ultimately, 
why does the structural phase transition appear to be 
fully developed in CuMn 02 , while in a-NaMn 02 it only 
manifests in a phase-separated state. Here, we an¬ 
swer these questions by combining numerical calcula¬ 
tions with local-probe experimental investigations. First, 
we determine the dominant intralayer exchange interac¬ 
tions by modelling the magnetic susceptibility via exact- 
diagonalization calculations. We demonstrate that the 
difference in the Tn for the two compounds can be under¬ 
stood via the mismatch of the two non-equivalent inter¬ 
chain exchange interactions, i.e., by the different extent 
of the frustration present in the two compounds, while 
the ME contribution is negligible in this respect. More¬ 
over, we provide the first experimental microscopic in¬ 
sight into the magnetism of CuMn 02 via nuclear 

magnetic resonance (NMR) and nuclear quadrupolar res¬ 
onance (NQR) measurements, as well as complementary 
muon spin relaxation (/iSR) measurements. These ex¬ 
periments clearly reveal that the ground state is more 
homogeneous in the Cu case than in the Na case and sug¬ 
gest that the stoichiometric CuMn 02 is on the verge of a 
phase-separation instability. Finally, our ab-initio calcu¬ 
lations suggest that the more homogeneous state of the 
stoichiometric CuMn 02 originates from an enhanced en¬ 
ergy difference (when compared to Q;-NaMn 02 ) between 
the two competing phases, born out of the magnetic- 
exchange and the elastic-energy changes below the Tn. 



T(K) T/J^ 


FIG. 1. Determination of J’s and Tn. (a) The mag¬ 
netic susceptibility y = M/H {M is the magnetization and 
H is the applied magnetic field) of CuMn02, measured by a 
SQUID magnetometer in a field of fioH — 0.1 T. The solid and 
dashed lines denote the best FTLM fits with the average and 
the scaled curves, respectively (see Methods for details). The 
former yields the exchange-coupling constants Ji = 53.5 K, 
J 2 /J 1 = 0.25, and the latter Ji = 52.1 K, J 2 /J 1 = 0.29. In¬ 
set shows the spatially anisotropic triangular spin lattice of 
the CuMn02 in the monoclinic setting, with intrachain Ji 
(thick bonds) and interchain J 2 (thin bonds) exchange con¬ 
stants. (b) The temperature dependence of the staggered 
susceptibility ys multiplied by Ji for spin-2 chains (adopted 
from Ref. l34f) . The solid line is a guide to the eye. The Neel 
transition temperatures Tn = 1.40Ji = 74 K in CuMn02 and 
Tn = 0.87Ji = 57 K in a-NaMn02 are predicted (dashed 
lines) by equation (p)l. 


RESULTS 

Determination of the dominant exchange interac¬ 
tions and the Tn. In order to understand the appar¬ 
ently significantly different properties of the two isostruc- 
tural compounds, CuMn 02 and Q:-NaMn 02 , a proper 
determination of the dominant terms in the Hamilto¬ 
nian is crucial. Therefore, we applied numerical finite- 
temperature Lanczos method (FTLM) simulations and 
density-functional theory (DFT) calculations. The for¬ 
mer were aimed at quantifying the two main magnetic 
exchange interactions (Ji and J 2 ) of the isotropic Heisen¬ 
berg model on the spatially anisotropic two-dimensional 
(2D) triangular lattice (inset in Figure [T^), 

= Ji^S,-8, +J2^Sfc-Sz. (1) 

(ij) [kl] 

Here, the first sum runs over the (stronger) intrachain 
bonds in one direction, while the second sum runs over 
the (weaker) interchain bonds in the other two directions 
on the triangular lattice of spins S = 2. 

The temperature-dependent magnetic susceptibility 
y(T) of this model was calculated for various m x n spin 


clusters (see Methods). However, even for the largest 
reachable cluster sizes (7 x 2), some finite-size effects re¬ 
main present at low T. With the use of results from many 
clusters {m = 2—7) these effects can, however, be reduced 
and the value of y at the thermodynamic limit is thus 
approached. Both, the average susceptibility curve from 
the two largest-size clusters (that gives a better approx¬ 
imation than the individual clusters; see Methods) and 
the susceptibility curve obtained in an approach similar 
to hnite-size scaling (see Methods) fit the experimental 
data above the Tn very well (Figure [TJi). A disagreement 
with the experimental data below the Tn, on the other 
hand, is expected, because the 2D Heisenberg model can¬ 
not account for a finite Tn. Both approaches yield very 
similar exchange-coupling constants, which we estimate 
to be Ji = 53 K and J 2 /Ji = 0.27(2). These are in 
very good agreement with recent ab-initio predictions^^, 
Ji = 56 K and J 2 /J 1 = 0.23, which, in principle, could 
be erroneous due to the unknown on-site repulsion^^. 

Our calculations thus confirm that the exchange inter¬ 
actions are indeed reduced for CuMn 02 compared to a- 
NaMn02, where^^ Ji = 65 K and J 2 /Ji = 0.44. Despite 
this fact, in CuMn 02 the Tn is increased with respect 
to that found in a-NaMn02 by more than 40%. So far. 
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this has been attributed to a difference in the interlayer 
coupling2^ J', which, however, is rather smal P^i^^ and 
can, therefore, only slightly affect the Tn on 2D Heisen¬ 
berg lattices^. Furthermore, the amount of frustration 
reflected in the J 2 /J 1 ratio should directly influence the 
Tn, as the frustration is known to suppress the spin 
correlations^^. In a-NaMn02 and CuMn02 the intra¬ 
chain exchange coupling is dominant. Therefore, these 
compounds can be regarded as systems of coupled spin 
chains, which is manifested in the one-dimensional char¬ 
acter of the magnetic excitations in Q;-NaMn02^. For 
such systems the Tn can be determined with the use 
of a random-phase approximatio n^^d*^ . In this approach 
the interchain coupling is treated at the mean-field level, 
whereas the intrachain interactions are treated exactly. 
For isotropic interchain coupling in a non-frustrated 
lattice the Tn is determined by the condition^&i^^— 
zJ 2 Xs{Tn) = 1 , where Xs{T) is the chain’s staggered 
susceptibility. Within this approach, we generalize the 
above condition for Tn to include the two interchain con¬ 
stants (J 2 a > J 2 h) pertinent to the triclinic phase of 
CuMn 02 and a-NaMn 02 , as well as the interlayer cou¬ 
pling J'; 

Xs(Tn) = 71—77- T N I , ( 2 ) 

k[z (J2a - J2h) + Z'J'I 

Here, z = z' = 2 corresponds to the number of 
neighbouring coupled chains and planes, respectively, 
while J 2 b adopts the minus sign because it frustrates 
the AFM order dictated by the larger J 2 a- The con¬ 
stant k renormalizes the coordination numbers and is re¬ 
duced from unit} d^d^ because of quantum effects^. As 
(d 2 a — J 2 h) I Ji = 0.083 and 0.035 in CuMn02 and a- 
NaMn02, respectively^^, while J'/Ji is expected to be 
an order of magnitude smalle r^^i^^ , we estimate the Tn 
by neglecting J' in equation ([2]) and by taking k = 0.7, 
which is appropriate for quasi-one dimensional cases (see 
Fig. 2 in Ref. [sl). This gives Tn = 74 K in CuMn02 and 
Tn = 57 K in a-NaMn02 (see Figure [Dd), which are in 
good agreement with the experimental values of 65 K 
and 45 K, respectively. We note that for anisotropic in¬ 
terchain couplings the constant k is expected to be fur¬ 
ther reduced and, ultimately, for— J' ^ 0 also fc —>■ 0, 
leading to Tn —>■ 0, which is consistent with the Mermin- 
Wagner theorem (no long-range order in the 2D Heisen¬ 
berg model at finite T). However, it has been shown^ 
that for quasi-2D systems the dependence of the Tn and, 
in turn also of k, on J' is sub-logarithmic. Therefore, for 
the exchange-coupling constants related to CuMn 02 and 
a-NaMn 02 , k will be somewhat, but not drastically, re¬ 
duced from the value 0.7, which is perfectly in line with 
the small theoretical overestimates of the Tn- 

This analysis reveals that the Neel transition is pre¬ 
dominantly determined by the Heisenberg Hamiltonian 
of equation O- Moreover, the ferrodistortive structural 
transition accompanying the magnetic ordering and lead¬ 
ing to the splitting of the two interchain exchange con¬ 
stants (J 2 a, J 2 b) in the triclinic phase is needed to ensure 



FIG. 2. DFT calculations, (a) The calculated difference 
in the total-energy density for the three different magnetoe- 
lastic components. The solid lines are linear fits for the xx 
and yy components, and a quadratic fit for the xy component, 
(b) The ab-initio calculated total energy of the relaxed mon¬ 
oclinic (m) and triclinic (t) structures of the CuMn 02 as a 
function of the volume for the antiferromagnetic (AFM) and 
non-magnetic (NM) cases. The insets zoom at the regions 
around the local minima of the relaxed structures. The global 
minimum of the energy is set to zero and the corresponding 
volume of the triclinic structure Vo is used for volume nor¬ 
malization. 


a finite Tn- The extent to which frustration is relieved 
in the triclinic phase of the CuMn 02 elevates its order¬ 
ing temperature above the ordering temperature in a- 
NaMn02. Other factors, such as the interlayer coupling, 
the magnetic anisotropy and the ME coupling, can, at 
best, only slightly shift the Tn- 

Total-energy change at the Tn. Having established 
that the magnetic ordering at the Tn is predominantly 
set by the 2D Heisenberg Hamiltonian and the tendency 
of both systems to remove magnetic degeneracy in the 
ground state by lattice deformation, the question that 
arises is what is the microscopic origin of such a complex 
transformation. In this respect, the ME coupling has 
been suggested as being the key facto r^^i^^i^^ in both a- 
NaMn02 and CuMn02. However, the ME coupling has 
been shown to be insubstantial in the former casein, and 
thus needs to be evaluated also in the CuMn02. The 
total-energy change at the Tn, associated solely with 
the magnetoelasticity, arises from the coupling terms^ 
bijeijTnimj between the strain-tensor components and 
the magnetization-direction vector m = {mx,my,mz)- 
The strength of the ME coupling, and consequently the 
corresponding contribution to the total energy, is propor¬ 
tional to the ME-coupling coefficients 6^ (i,j = x,y,z). 
The coefficients bxx = 2.3 MJ/m^, byy = 1.6 MJ/m^ 
and bxy = 3.4 MJ/m^ are determined as linear terms 
in the calculated dependence of the total-energy-density 
change Afij on the strain (see Methods), which is shown 
















4 



FIG. 3. NMR results, (a) 80-K NQR and NMR (inset) spectra of CuMn02. The solid lines represent a simul¬ 
taneous NMR/NQR fit (see Methods for details), assuming a Gaussian distribution of the NQR frequencies = 

-I- r;2/3, that yields the quadrupolar frequency = 27.0(1) MHz, the asymmetry parameter rj = 0.20(5), the 
isotropic hf shift Kht = 1.6(1)% that is much larger than the dipolar shift Kd — 0.11%, and the individual line widths 
= 0.17(1) MHz, = 0.14(1) MHz. The dashed lines show the center of the NQR lines and the reference NMR frequencies 
corresponding to a zero magnetic shift. The ®®Cu NQR spectrum at 4.6 K is added for comparison, (b) The temperature 
dependence of the ®®Cu NQR line width ®®5 and the line position ®®r'c. The inset highlights the hf paths through the 
sites that provide the coupling of each Cu nuclei with six surrounding Mn®^ magnetic moments (arrows), ordered with the 
magnetic wave vectop^i k = (0, ^). (c) Gomparison of the temperature-dependent ®®Cu NQR/NMR spin-lattice relaxation 

rate l/Ti in the CuMn02 and ^®Na NMR l/Ti in the Q:-NaMn02 (Ref. 0 ). The latter is normalized by the squared ratio of 
the hf coupling constants. The error bars represent the standard deviation of the fit parameters. 


in Figure [2^. The ME energy gain is the largest for the 
shear-strain component exy, which is associated with the 
monoclinic-to-triclinic deformation. However, for the ex¬ 
perimental strain e^y = 0.0028 (see Methods) the mag- 
netoelastic energy change amounts to only 2.7 peV per 
triclinic unit cell, which is very similar to the value of 
2.5 peV found in a-NaMn02. 

The contribution of the ME coupling to the total en¬ 
ergy change at the Tn is thus negligible in both com¬ 
pounds. Therefore, the complex phase transition at the 
Tn has to reflect changes in the magnetic-exchange and 
elastic energies^^. Our ab-initio calculations of both rele¬ 
vant contributions to the total energy in the CuMn 02 for 
the non-spin-polarized case, relevant to non-magnetically 
ordered structures, reveal that the monoclinic structure is 
energetically lower than the triclinic one, although only 
by ~1 meV per 4 formula units (f.u.); see Eigure 
This is in-line with the (72/771 crystal symmetry found 
experimentally at room temperature. However, once the 
magnetic order sets in, the total energy of the triclinic 
structure is lowered below that of the monoclinic struc¬ 
ture by about 4 meV per 4 f.u. This change of ^1.25 meV 
per Mn®+ ion is mainly a consequence of the exchange- 
energy decrease during the structural phase transition, 
associated with the removal of the degenerate magnetic 
states due to the interchain frustration. The result¬ 
ing splitting of the interchain exchange constants by2^ 
'ha. — 72b = 0.4 meV releases S^{J 2 a — 72b) = 1-6 meV of 
energy per Mn that is slightly larger than the total-energy 


change at the Tn, as it is partially spent to compensate 
for the elastic-energy increase in the triclinic phase. 

The calculated total-energy difference below the Tn of 
1 meV per f.u. in CuMn02 between the two structures 
that is about 3-times above the calculation error bar (see 
Methods), is markedly larger than in the Q!-NaMn02, 
where the calculated difference was below the calcula¬ 
tion error bar; <0.5 meV per f.u4^. However, in abso¬ 
lute terms this difference is small, even in the CuMn 02 , 
so that a competition between the near-degenerate mon¬ 
oclinic and triclinic structures is expected for both com¬ 
pounds. Experimental local-probe magnetic techniques 
are then essential for highlighting possible differences be¬ 
tween the two systems. 

NMR/NQR insight to the magnetism. Informa¬ 
tion about the magnetic properties of CuMn 02 on the 
local scale are revealed in the NQR/NMR experiments 
via the hyperfine (hf) coupling Hhf of the electronic 
and the nuclear magnetic moments. Moreover, 

the quadrupolar splitting in the electric-field gradient 
(EFG) provides information about the material’s struc¬ 
tural properties. The NQR spectra measured in zero 
field correspond to a single line for each copper isotope^, 
while the powder NMR spectra are structured (see Meth¬ 
ods for details). Our simultaneous fit of the NQR 
and NMR data at 80 K (Figure [5^) yields a hf cou¬ 
pling constant = 2.3(1) T/ps that is significantly 

larger than the coupling constant = 0.11(1) Tjfis 
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foundi^ in Q:-NaMn02. Since Ahf scales with the or¬ 
bital overlap, the charge transfer from the Mn^+ ions 
to the interlayer cations (Cu+ or Na'*') is much larger 
in the CuMn02 than in the Q:-NaMn02. This implies a 
stronger J' in the former compound and thus is in line 
with the somewhat better agreement between the experi¬ 
mental and the predicted Tn, as the renormalization fac¬ 
tor k in equation ([2|) is closer to the used value of 0.7. 

The width 6 of the NQR spectra at 80 K, amount¬ 
ing to 0.16 MHz/27.2 MHz = 0.6% of the line-position 
value r'c already in the paramagnetic phase (Figure [Sb), 
is rather large. The line widths reveal spec¬ 

tral broadening being in accordance with the quadrupole 
moments and contradicting the gyromag- 

netic ratios Therefore, sizeable structural 

distortions of the local environments must be present. 
The temperature dependence of both Vc and 6 shows a 
pronounced sudden increase below the Tn (Figure [SJd), 
clearly marking the phase transition. The anomaly in Vc 
at the Tn is attributed to the structural transformation 
of the CuMn 02 sample, directly affecting the quadrupo- 
lar frequency z/q. Namely, static internal magnetic helds 
below the Tn cause a symmetric broadening/splitting of 
the NQR line so that its center of gravity is unaffected^. 
On the other hand, the pronounced increase of <5 by a fac¬ 
tor of ^2 at the Tn, exceeding the change of i/c by several 
orders of magnitude, can only be magnetic in origin. 

We must emphasize that the existence of the NQR 
signal below the Tn is unexpected. Namely, in the 
frame of the homogeneously ordered magnetic phase^l 
with k = (0,^,^) the Cu nuclei would experience ex¬ 
tremely large internal magnetic fields. Although the Cu 
site is a structural center of inversion, the magnetic or¬ 
der breaks this symmetry, as the spins at ±r from a 
given Cu site are FM ordered (see the inset in Fig¬ 
ure ISb), in contrast to the a-NaMn02, where the or¬ 
der of the two corresponding spins is AFMi^. Such a 
spin configuration in CuMn 02 yields a large local hf field 
Hhf = ®^Ahf/3 X /r ~ 2.2 T (/i = 3.05/iB is the size of the 
ordered^^ Mn^+ moment). This field leads to extremely 
broad NQR spectra, = ®^7/27r x 2Hhf — 50 MHz, 
being two orders of magnitude broader than the experi¬ 
mental ones. Indeed, the NQR signal below the Tn corre¬ 
sponds to a minority fraction of all the ®^Cu nuclei, while 
a majority of the signal is lost at the Tn due to the onset 
of large internal fields. Namely, the Boltzman-corrected 
intensity of the NQR signal at 4.6 K, when further cor¬ 
rected for nuclear relaxation effects, is smaller than the 
intensity at 80 K by a factor of ^ 17. This reveals that, 
unexpectedly, about 6% of all the Cu sites in our sam¬ 
ple experience small or no internal magnetic fields and 
do not correspond to the reported homogeneous mag¬ 
netic phase. We note that the AFM order of the mo¬ 
ments positioned symmetrically with respect to the Cu 
site, or the absence of any order, result in a zero static lo¬ 
cal magnetic field at the Cu site, and would explain the 
NQR-observable sites below the Tn. Since this minor¬ 


ity signal exhibits clear anomalies at the Tn (Figure [Sh) 
it is obviously well coupled to the bulk that undergoes 
the magnetostructural transition. This is confirmed by 
the temperature dependence of the spin-lattice relaxation 
rate, 1/Ti, that shows a maximum at the Tn due to criti¬ 
cal spin fluctuations^ related to the magnetic instability 
of the bulk. However, in contrast to the monotonic de¬ 
crease found in the Na-based compound below the Tn, in 
the CuMn02, another clear maximum in both the NMR 
and NQR 1/Ti is observed at around 10 K. This reveals 
an, as yet, unobserved instability that could be either 
magnetic or structural in its nature. 


Probing the magnetic disorder with p,SR. In or¬ 
der to provide more insight into the magnetic state in 
the CuMn 02 that is, according to the unexpected mi¬ 
nority NQR signal, apparently not as homogeneous as 
inferred from previous bulk measurements, we resorted 
to the ^SR local-probe technique. Moreover, this tech¬ 
nique reveals details about the low-temperature anomaly 
in the NMR/NQR relaxation at 10 K. In contrast to 
NQR/NMR, which is limited because the intrinsic signal 
disappears below the Tn, the /rSR measurements can as¬ 
sess the magnetic properties of the entire CuMn 02 sam¬ 
ple also below the Tn. This time a hf/dipolar coupling 
between the electronic magnetic moments and the muon 
magnetic moment is utilized after a muon stops in the 
sample. The resulting local magnetic field at the 
muon site affects the ^SR asymmetry A{t) that is pro¬ 
portional to the muon polarization processing in B^. 

In CuMn02, the weak-transverse-field (wTF) exper¬ 
iment that effectively keeps track of the temperature- 
dependent ordered part of the sample by measuring the 
amplitude of the oscillating ^SR signal^, reveals that 
the fraction of the muons detecting large frozen internal 
fields starts growing already below 80 K (Figure 0^); i.e., 
far above the Tn, which can be attributed to develop¬ 
ing short-range spin correlations^. Below the Tn, the 
whole sample becomes magnetically ordered within only 
a few kelvins, leaving no room for a non-frozen fraction 
above the experimental error bar of a few percent. Simi¬ 
lar information is obtained from the zero-field (ZF) /iSR, 
where the initial asymmetry strongly decreases around 
the Tn and at low temperatures reaches 1/3 of its high- 
temperature value (inset in Figure|4]3). Such a reduction 
is characteristic of the establishment of strong static in¬ 
ternal fields in powder samples. Statistically, in 1/3 of all 
cases the muon magnetic moment is aligned parallel to 
the B^ and therefore exhibits no precession, while rapid 
oscillations of the asymmetry in other cases diminish the 
/iSR signal on a coarse time scale. 

A detailed look at the ZF relaxation curves below Tn 
(F igure 01:) also allows for the detection of the quickly- 
oscillating component. Similar to the a-NaMn02 casein, 
the two-component model 
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FIG. 4. fiSH. results, (a) The temperature-dependent magnetically ordered volume fraction [1 — y4o(r)/Ao(120 K)] of 
the CuMn 02 , derived from the wTF /rSR asymmetry Awtf data (inset); solid lines are fits to the model AwTF(t) = 
^o(T) cos (y^BwTpt) -I- C{T), where 7 ^ = 27r x 135.5 MHz/T is the muon gyromagnetic ratio, BwTF is the transverse 

applied magnetic field and At the transverse muon rela:xation rate. The Ao{T) term corresponds to muons experiencing no 
sizeable static internal magnetic field, while the C{T) term describes those muons that reside at sites with large static fields 
3> Bwtf). (b) The longitudinal muon relaxation rate derived from the ZF muon asymmetry Azf data (inset); solid lines 
are fits to the model AzF{t) = A'o{T)e~^^^, where the initial asymmetry Ao{T) is temperature dependent to account for the 
disappearance of the oscillating component below the Tn. (c) /rSR asymmetry of CuMn 02 (upper panel) and Q:-NaMn 02 (lower 
panel; adopted from Ref. 13 ) at 5 K. The solid and dashed lines represent the corresponding fits to the "two-component" model 
of equation m and a model with only one oscillating component, respectively. Fitting to the CuMn 02 data yields = 0.71 and 
1.71 for the former and the latter models, respectively. The error bars represent the standard deviation of the fit parameters. 
For the muon asymmetry data the latter are set by the square root of the total number of detected positrons. 


fits well with the experimental data. Here, fj denotes 
temperature-independent probabilities that the muons 
stop at either of the two magnetically non-equivalent 
stopping sites j. The preferential site is occupied in 
70(5)% cases. The internal field at this site is only 
slightly higher than at the second site (0.59 and 0.54 T 
at the first and the second site, respectively, at 5 K); 
however, a fit with only a single oscillating component 
(the dashed line in Figure |4}:) results in a much poorer 
agreement with the data. The damping rate of the oscil¬ 
lations Xtj that is due to the finite width of the local-field 
distributionsi^ in CuMn 02 is reduced by a factor of ^3 
when compared to the a-NaMn02 (Figure[4j:), indicating 
more homogeneous magnetism. 


In the ZF experiment, the magnetic phase transition 
at the Tn is expressed as a maximum of the longitudi¬ 
nal muon-relaxation rate Al, like in the NQR/NMR re¬ 
laxation experiments. Moreover, the second maximum 
observed in the NQR/NMR experiments at 10 K is also 
found in the /rSR. Since the ZF pSR signal corresponds 
to the total volume of the sample and the muons are 
only sensitive to magnetism, this reveals that the low- 
temperature anomaly is of magnetic origin and is intrin¬ 
sic to the CuMn 02 system. 


DISCUSSION 

The FTLM and DFT numerical calculations provide 
a solid basis for addressing the experimentally observed 
similarities and differences between the CuMn 02 and the 
a-NaMn02. Considering the latter calculations in the 
magnetically ordered state, the triclinic phase is energeti¬ 
cally preferred in both compounds. With increasing tem¬ 
perature, the staggered susceptibility decreases and this 
leads to a finite Tn. Above the Tn the exchange-energy 
gain associated with the magnetically ordered state dis¬ 
appears, which in turn leads to a structural transforma¬ 
tion to the monoclinic phase that is energetically pre¬ 
ferred in the non-magnetic state. The isotropic Heisen¬ 
berg Hamiltonian of the spatially anisotropic triangular 
lattice is dominantly responsible for elevating the Tn in 
the CuMn 02 with respect to the Q;-NaMn 02 , while the 
magnetoelastic and the interlayer couplings play a less 
important role. 

On the other hand, our local-probe experiments on 
the CuMn02revealed some subtle, yet profound, fea¬ 
tures that should be carefully considered in the attempt 
to understand the presence/absence of nano-scale phase 
separation in the spatially anisotropic triangular lat¬ 
tice. Both, the NQR/NMR and the pSR investiga¬ 
tions demonstrated that CuMn 02 undergoes a magne- 
tostructural phase transition at Tn = 65 K almost com¬ 
pletely. The minority NQR component (^6%) that re- 
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mains present below the Tn can be explained by re¬ 
gions where the interlayer magnetic ordering is FM in¬ 
stead of being AFM, as the latter causes the disappear¬ 
ance of the NQR signal due to large local fields. The 
NQR signal exhibits a magnetic anomaly around 10 K, 
which is expressed by the increased relaxation rates of the 
NQR/NMR as well as the fj,SR. Since /rSR, on the other 
hand, detects a bulk magnetic signal, the small NQR 
component is apparently coupled to the bulk magnetic 
phase. This is further confirmed by the line position and 
the width of the NQR spectra, changing considerably at 
the Tn. The coupling with the bulk phase can then be 
regarded in the context of the nano-scale phase inhomo¬ 
geneity. A comparison of the ZF /iSR asymmetry curves 
of the CuMn 02 and the a-NaMn 02 is quite informative 
in this respect. The notably reduced damping of the os¬ 
cillations in the ordered phase of the former compound 
provides evidence of much narrower field distributions, 
and hence less disorder. This conclusion is also in line 
with the number of the interlayer cation (Cu+ and Na+) 
sites experiencing internal fields that do not comply with 
the symmetry of the bulk magnetic order, which in the 
CuMn02 is decreased to 6%, from the 30% foundi^ in 
the a-NaMn 02 . 

The magnetostructurally inhomogeneous ground state 
of the a-NaMn 02 on the nano-scale has previously been 
attributed to the combined effects of geometrical frustra¬ 
tion and near-degenerate monoclinic and triclinic struc¬ 
tural phasesi^. We believe that the key factor control¬ 
ling such an inhomogeneity is the difference in the to¬ 
tal energy of the two competing phases in the magnet¬ 
ically ordered state. This difference is notably larger 
in the CuMn02 (1 meV per f.u.) than in the a- 
NaMn 02 , where it is below the computational error bar 
(<0.5 meV per f.u4^). In the latter compound, an in¬ 
finitesimal quenched disorder, locally favouring one phase 
over the other, can then be held responsible for trigger¬ 
ing the phase separation. Similar effects are suppressed 
in the stoichiometric CuMn 02 , but would become en¬ 
hanced for larger deviations from perfect system uni¬ 
formity. Indeed, enhanced strain, acting as a precursor 
of the monoclinic-to-triclinic structural phase transition, 
has been observedi^^ in the high-temperature mono¬ 
clinic phase in stoichiometric Q;-NaMn 02 , while in the 
CuMn02 a Cu-Mn off-stoichiometry is required to pro¬ 
duce such a strain2^. Moreover, the diffuse magnetic scat¬ 
tering characteristic of 2D correlated regions that coexist 
with sharp magnetic Bragg peaks (one of the signatures 
of the inhomogeneityi^ found in the Q!-NaMn 02 ) is also 
found^^i^OiSi in the CuMn02. However, in contrast to 
the a-NaMn 02 , where it persists to low temperatures, 
in stoichiometric CuMn 02 it gradually gives way to the 
3D ordered phase below the Tn. Interestingly though, 
in off-stoichiometric samples^ the volume fraction of the 
2D-correlated phase shows no decrease below the Tn, im¬ 
plying that the 2D-ordered regions keep competing with 
the 3D order at low temperatures. The total-energy dif¬ 
ference of the competing phases below the magnetostruc- 


tural transition, reflecting the interplay of the magnetic- 
exchange and the elastic energies, then seems to deter¬ 
mine the amount of disorder required to stabilize the in¬ 
homogeneous ground state on a geometrically frustrated 
triangular lattice. Systems with near-degenerate compet¬ 
ing phases can be locally perturbed more easily. Such an 
inhomogeneity may, therefore, be a more general feature 
of geometrically frustrated magnets. 


METHODS 

Finite-temperature Lanczos method simulations. 

Calculations of the spin susceptibility for the S = 
2 Heisenberg model on the anisotropic triangular lat¬ 
tice (equation H])), were performed with the finite- 
temperature Lanczos method (FTLM) "*^^i^° and were used 
to determine the leading exchange couplings Ji and J 2 
in the CuMn02. Within the FTLM finite-size clusters 
are diagonalized in a similar manner as for the standard 
exact diagonalization Lanczos method (at T = 0) and 
additional random vector averaging over the R vectors 
is employed to determine the properties at T > 0. Typ¬ 
ically, i? ~ 10 suffices for the largest systems and the 
lowest T, while smaller systems require a larger R. The 
limitations of the method are mainly set by finite-size 
effects, which are the largest at low T and determine 
the lowest reachable T(^ Ji). In order to reduce the 
finite-size effects we used periodic boundary conditions, 
adjusted cluster shapes, the largest reachable cluster sizes 
(up to = 14 sites), and additional approximations for 
the values in the thermodynamic limit. 

The temperature-dependent magnetic susceptibility 
x(T) of the model given by equation ([T|) was calculated 
previously in Ref. for Q;-NaMn02. It was shown that 
in the regime of interest, elongated spin clusters are the 
most appropriate. In particular, if a cluster has m in¬ 
dependent spins in the Ji direction and n spins in the 
J 2 directions, it was realized that due to J 2 < Ji and 
two competing J 2 bonds, y does not depend on n for 
n > 2 (see Fig. 2 in Ref. [3^. This fact allows us to 
reduce the finite-size effects by using a larger m. We 
note that due to the alternating behaviour of x with 
m (see Figure [5]), which originates in periodic boundary 
conditions and antiferromagnetic spin-spin correlations, 
the average susceptibility curve from the two largest-size 
clusters (m = 6 and 7) is a better approximation than 
the m = 7 curve. 

Scaling-like approximation for the susceptibility. 

The results for several different sizes of finite clusters and 
their systematics, shown in Figure [51 allow a scaling-like 
analysis to obtain a better approximation of y in the 
thermodynamic limit. Typical scaling analyses use scal¬ 
ing functions of the form x(A) = a + b/N and additional 
higher terms when needed; e.g., c/N"^. Since our cal¬ 
culations are limited to rather small maximum system 
sizes by S' = 2, we also use the results from small sys- 




FIG. 5. The flnite-size effects in the FTLM calculation. 

The calculated susceptibility on various mxn clusters for the 
optimal parameters Ji = 53 K and J 2 IJ 1 ~ 0.27. The curve 
averaging the 6x2 and 7x2 data and the scaled curve are 
also shown. 


terns (starting with m = 2). Consequently, such scaling 
functions are not appropriate. In particular, at higher T 
(see, for example, T > 5 Ji in Figure [S]) x has already 
converged with N for systems with m > 5, while for 
m < 5 notable finite-size effects are seen. Therefore, the 
scaling function should be close to a constant for 1/iV 
smaller than some value, while at larger 1/N, the scal¬ 
ing function should allow for a stronger N dependence. 
For these reasons we use a generalized scaling function 
of the form x(iV) = a -I- 6[exp(c/A^) — 1], which corre¬ 
sponds to typically used functions in the limit of small 
1/A^. In order to also capture the alternating compo¬ 
nent of X with m (Figure [S]) we add, in a similar fash¬ 
ion, the term bi [exp(ci/A^) — 1] (—1)^/^. Such a scaling 
function also gives a correct (converged with N) result 
for high T, while typical scaling functions fail in this re¬ 
spect. We have performed such a scaling for each T sep¬ 
arately. However, since we are limited to small systems 
with notable finite-size effects at low T (< di) and since 
the scaling function has many parameters, the result of 
such an analysis should not be taken as a strict thermo¬ 
dynamic limit. Rather, it should be regarded as a next 
approximation of it, compared to the result from simpler 
averaging of the two largest-cluster curves. 

Density-functional theory calculations. The cal¬ 
culations of the total energies and the magneto-elastic 
(ME) coupling coefficients were performed within the 
framework of the density-functional theory (DFT) and 
the generalized-gradient approximation (GGA)^ for the 
exchange-correlation contribution by applying the Quan¬ 
tum Espresso code^. The electron-ion interactions were 
described by the Vanderbilt ultrasoft potentials^ includ¬ 
ing the spin-orbit coupling for the Mn atoms. The plane- 
wave cut-off parameters were set to 585 eV and 4678 eV 


for the expansion of the wave functions and the potential, 
respectively. In order to take into account the proper an¬ 
tiferromagnetic ordering, the 1x2x2 supercells of the 
monoclinic and the triclinic structures were used. The 
calculations of the total energies as a function of the 
unit-cell volume for the different types of magnetic or¬ 
dering were carried out by using 4x8x2 reciprocal 
vectors in the full Brillouin zone (BZ) for the Methfessel- 
Paxton sampling^ integration. The criterion for the self 
consistency was the total-energy difference between two 
subsequent iterations being less than 10“®Ry. The mon¬ 
oclinic phase was further optimized by minimizing the 
total energy and the inter-atomic forces with respect to 
the lattice parameters and the atomic positions. The re¬ 
sulting structure served as the zero-strain reference for 
the calculations of the ME coefficients that are based 
on the evaluation of the total-energy differences of the 
order of < 10“"^ Ry, which is also the accuracy for the 
determination of the total-energy differences between the 
monoclinic and triclinic phases in Figure calculated 
per 4 f.u. The tests yielded 8x16x4 reciprocal vectors 
in the full BZ to be enough for well-converged results. 

Determination of the magnetoelastic coupling. 

The magnetoelastic coupling constants bij are calculated 
from the associated magnetoelastic energy density. This 
contains the products bijeijrriimj of the strain-tensor 
components and the components rrii of the normalized 
magnetization. The form of the magnetoelastic energy 
density is determined by the symmetry of a particular 
system^. The magnetism of the GuMn02 and the a- 
NaMn 02 is essentially two-dimensional; therefore, only 
the terms with the lateral strain-tensor components are 
important. For the monoclinic symmetry, these include 

^xx^xx'^x^ 
bxy^xy^x'^y 5 

byyeyymy. (4) 

Individual magnetoelastic terms are then determined by 
specifically choosing strain components and magnetiza¬ 
tion directions and calculating the total-energy density 

f (^^xx 5 ^yy 5 ^xy ; '^x : '^y 7 '^z ) : from 

Afxx = f(eo, 0,0,1,0,0) - f(eo, 0,0,0,0,1), 

A/,, = /(O, eo, 0,0,1,0) - /(O, eo, 0,0,0,1), 

Nxy = /(0,0, eo, 1,1,0) - /(0,0, eo, 0,0,1). (5) 

The above total-energy differences are calculated ab- 
initio as a function of eo for relaxed crystal structures. 
In GuMn02, Afxx and Afyy change linearly with in¬ 
creasing strain at least up to eo = 0.03, while an addi¬ 
tional quadratic term is observed in Afxy (Figure [21i). 
The experimental strain value Cxy = 0.0028 that is ob¬ 
tained by calculating the relative shift of the Mn^'*' ions 
in the triclinic structure, when compared to the mon¬ 
oclinic structure (based on high-resolution synchrotron 
XRD data analysis!^), is an order of magnitude lower. 
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Therefore, the linear term is dominant for all three con¬ 
tributions and allows the extraction of the three magne¬ 
toelastic constants b^x = 2.3 MJ/m^, byy = 1.6 MJ/m^ 
and bxy = 3.4 MJ/m^. 

Nuclear magnetic/quadrupolar resonance. 

(/ = 3/2) NMR/NQR measurements were performed 
on a high-quality powder sample with the same phase 
purity and stoichiometry as in the study presented in 
Ref. [^. The NMR/NQR spectra and the spin-lattice 
relaxation were measured between 4.6 K and 120 K in 
a magnetic field of 8.9 T (NMR) and in zero magnetic 
field (NQR) on a custom-built spectrometer. Frequency 
sweeping and a solid-echo pulse sequence were used for 
recording the spectra, while a saturation recovery method 
was used for measuring the spin-lattice relaxation. Typ¬ 
ical 7 r/ 2 -pulse lengths were 3.5 ps and 6 ps in the NMR 
and NQR experiments, respectively. The reference NMR 
Larmor frequencies of = 100.728 MHz and = 
107.908 MHz were determined with a O.IM NaCl-solution 
reference by taking into account the gyromagnetic ratios 
237 = 27r X 11.261 MHz/T, ^37 = 2 tt x 11.295 MHz/T 
and 6^7 = 27r X 12.089 MHz/T. 

The NQR spectrum of each isotope is particularly sim¬ 
ple, as it is given by a single line^ at = 

63,65j^q^ 1 _j_ with the ratio of the quadrupolar fre¬ 
quencies ^^vci/^^vci = 1.08 fixed by the corresponding 
quadrupolar moments and the EFG tensor asymme¬ 
try parameter being 77 = {Vxx — V/y)/t4z. The NMR 
spectrum is more complicated, because the applied mag¬ 
netic field Bq breaks the symmetry in the spin space. 


The central-transition (—1/2 o 1/2) powder NMR line 
adopts a characteristic structure because of the angular- 
dependent NMR shift K from the reference frequencies 
uo = jBo/2tt, 63.65;^ = = Km + iFd + 

In analogy^ to the a-NaMn02, we take the hf shift 
Km = is the reduced Planck constant) to 

be isotropic, while the dipolar contribution Tfd and the 
quadrupolar— shift can be accurately calculated. 

The former has a uniaxial symmetry and is calculatecU^ 
(iFd = 0.11% is the dominant eigenvalue) by taking into 
consideration all the Mn3+ paramagnetic spins around 
a given Cu site within a sphere large enough to ensure 
convergence. 

A homogeneous life-time broadening of the NQR spec¬ 
tra is negligible. The spin-spin relaxation time = 
46 /is at 80 K yields = 6.9 kHz, which is much 
smaller than the spectral width. The spin-lattice relax¬ 
ation is of magnetic origin. We find the isotopic effect 
®3 Ti/®3Ti= 0.86 that is in accordance with magnetic re¬ 
laxation dictating 1/Ti oc 7 ^A^£. 

Muon spin relaxation. The /iSR investigation was 
carried out on the General Purpose Surface muon (GPS) 
instrument at the Paul Scherrer Institute, Villigen, using 
the same powder sample as in the NMR/NQR experi¬ 
ments. Zero-field (ZF) and weak-transverse-field (wTF) 
measurements in a 3 mT magnetic field were performed 
in the temperature range between 5 and 120 K. The veto 
mode was utilized to minimize the background signal. 
The ZF pSR measurements below the Tn revealed that 
each muon stops at one of the two possible non-equivalent 
stopping sites, like was observecU^ in a-NaMn 02 . 
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